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/- VOLUME DISCHARGE IN A GAS, EXCITED BY AN ELECTRON BEAM '/ 357 * 

/' -under CONDITIONS OP NONUNIPORM IONIZATION 


0. B. Evdokimov-, G., A. Mesiata, and V. B. Ponomarev** 

I. INTRODUCTION : 

Volume discharge in gases excited by electron beams [1] has found 
wide applications in powerful lasers [2] in powerful commutators [3],’ 
plasmatrons, etc. The theory of such discharges assumes that the dis- 
tribution in the rate of generation of the electron-ion pairs is uni- 
form along the depth of the discharge gap. However, in reality this is 
not always the case. Because of the probabilistic nature of electron 
scattering in gases, the spread in the energy spectrum of the electrons 
upon their passage through the foil and into the gas , and also because . .. 
of time-variations in. the voltage of the accelerating diode, the dis- 
tribution of , the generation rate of the electron-ion pairs can turn out 
to be non-uniform. In effect, the beam electrons become thermalized 
[4]. This leads to nonuniform concentration of plasma in the discharge , 
column and to a nonuniform electric field, which in certain cases leads ' 
to instabilities in volume discharge [53. Since these instabilities 
are connected with the special properties of the Injected beam, they 
are called injection instabilities [53^ In the case of lasers, the non- 
uniform ionization leads, in addition, to a nonuniform population in- 
version which is manifested in the optical steady-state properties of 
the laser. The question of distribution of the electric field under 

*Numbers in margin indicate pagination .in jforeign text . 
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the conditions of nonuniform ionization is discussed in part and on a 
qualitative level in [6]. In contrast with [6] we will derive the 
steady-state electric field distribution taking into account shock 
ionization. We will also consider the voltage-current and energy cha- 
racteristics and will take up the question of how the discharge pro'p-' 
erties, at short-time intervals, of injected beams are affected by 
the nonuniformity in the rate. of degeneration of the electron-ion 
pairs . 

We will start with the following system of equations for the one- 
dimensional case.: 

^'^-^Z =\y+ctn-v.—^n-n++g, 

I dt dx 

— 1 + — t_L=^+ort-y--pn-n+, CD 

dtdx 

dE e . , r 

ax e ■ 

Here, 'i' is the volume rate of generation of electron-ion pairs, the 
niAmber of electron-ion pairs created by the fast electron beam by 
ionization . of atoms per second per. unit volume; g is the volume speed . 
of /injection of beam electrons as a result of thermalization; n^, u^, 
y^'are the concentration, mobility, and . speed of the ions and electrons 
/respectively; e is. the electron charge; E. the electric field Intensity; 

'3 - recombination coefficient; a the coefficient of shock ioniza- 
tion. Here and in what follows, the anode is located at x = 0 and the 
cathode - at x = d. . • ■ ,/358 

By order of magnitude 'g/y ’^ei/(,Ty , where is the average energy of 
electron-ion pairs created, <f> is the average energy of fast electrons. 

We assume that the thermalization time of fast electrons is short. 


II. Electric Field Distribution in the. Steady State . 

. One of the remarkable properties of volume discharge at large ’}' is 
the small value of the cathode voltage drop Uj^, in comparison with the 
applied voltage U , U, << U [7]. Let us consider a discharge column 

O Jx -O 

where we neglect ion drift. In this case, our system of equations (1). 

2 ' ./■ ■ : 



■ . becomes : 

W (z)+g(x)+aii-n-E=^n+n-, 
dja ■ 

/«=/o(*)+en-n-E, ~r~‘=eg{x), 

ax ‘ 

dE e(ra+— n_) 

dx 8 ’ 


( 2 ) 

(3) 

(4) 


where j is the beam current density through the gas, at a depth X- In 
the steady-state case, usually j^(x) << ey_n_E; hence, j ey_n_E. Let 
us suppose, just as in the case of uniform ionization, that in the dis- 
charge column n;^ is approximately equal to n_, i.e.,, ■ 


!«+— n-|<n-, 

.Then, by eliminating the electron concentration from (3), we obtain in 
an. implicity form. the field distribution 


e'V (x) +«(£)/ “ » 


( 6 ) 


where we have set y_ = ,y. Figure 1 gives examples of typical distri- 
butions 'i'(x) obtained by solving the kinetic equation which describes 
the transport of monoenergetic electrons in nitrogen in a uniform field 
and by taking into consideration the separating foil. The fast elec- 
trons are Injected through the cathode. The function 'i'(x) decreases 
from the. cathode to the anode; nowhere is it equal to zero, since an 
electron beam propagating through the medium is accompanied by brehms- 
strahlung. Approximately ¥ _ 'v 10 t 10 . With increasing 

ITlXri ITl&X 

field intensity y(x) becomes more uniform. 


Clearly, starting ,wlth some minimum value E(x = xm) = near 

the cathode, the- field .increases in the direction of the anode. If 
dW 

— — >0 for all X, then x^ is at the boundary between the discharge co- 
lumn and the cathode region. When the discharge voltage is several 
times smaller than the static breakdown voltage, then at the point 
where the electric field is a minimum, the inequality YCx^) >>• o(E^^^)j/e 
holds and therefore 






p/* 




(7) 
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and by combining (6) and (7) for 0 < x < we have: 

, / . y(x) _ E min 4Lgmi..ce(g) 


( 8 ) 


It is possible to show, by considering (8) and (4), that condition /359 
(5) is true if 

(9) 

2cH'‘(x) dx \ ‘ 


Condition (9) takes into account the fact that E<E^,JW [x„)/'^V (x) . 

Hence, whenever shock ionization is significant, condition (9) is an 
overestimate and is fulfilled for realistic values of Y(x) (see Figure 
1). Figure 2 shows an example of. a calculation of the field E(x) for 
(ilscharge in nitrogen. It shows the influence of shock ionization on 
'the distribution of the field. The curve labeled 2 is a calculation 
without shock ionization. Curve 3 is a calculation taking shock ioni- 
zation into account according to equation (8), with parameters y = 

500 cm^/sec; 3 = 2 k lO”*^, the. a(E) dependence is taken from [8]. 




Figure 1: The distribution of the Figure 2: The distribution of the 

generation speed of electron-ion ■ ■ electron-ion pair generation rate 
pairs in nitrogen at atmospheric (curve 1) in relative units in 
pressure behind an aluminum foil ; nitrogen at atmospheric pressure 
of 50 y thickness. Curve 1 - in- at beam current 5*10”^a, initial 

put energy 150 keV, E/p 'v 4 kV/ energy 200 keV, foil thickness 

atm* cm; Curve 2 - energy 200 keV 50 y; 2-^ electric field distri- 
E/p = 2 kV/atm*cm; Curve, 3 - butlon without taking into ac- 

energy 200 keV, E/p = 4 kV/atm*cm count shock ionization; 3- elec- 

at the beam point 5. * 10”^ a/cm^. trie field distribution with shock 

ionization. 


The calculation uses the distribution Y(x) (curve i) calculated 

"min* - /^min^ 


for a uniform field E - E„^„. Since E(x) ^ the real distribution 
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is a smoother function and hence, because of Equation (8), curves 2 
and 3 give an upper limit for the field distribution. 


3. CURRENT- VOLTAGE AND ENERGY CHARACTERISTICS 

The relationships (6), (7), and (8) allow us to calculate the cur- 
rent voltage (CVC) and the energy characteristic of volume ^as dls-. 

t 

charge with Inhomogeneous Ionization with the condition j 

where Is the applied gap voltage. If We neglect shock-lonlzatlon,: 

then f^rom_(6) It follows that 

- ■ - — 


Eix)^ 


Integrating (10) from x 


0 to x^, 
m 


we obtain CVC 


( 10 ) 


7360 




enUt 


( 11 ) 




here. 




-J 

X, •' 




dx 


( 12 ) 


V'F(x) 


• When U^/d Is Increased, the distribution 'i'(x) far away from the 
cathode becomes more uniform (Figure 1) and according to (11) and (12) 
for sufficiently large distances between the electrodes, we can expect 
that with Increasing U , the current will increase more rapidly than 
linearly. 



JO d,CM 


' , Figure 3:. Distribution of genera- 
tion rate 'l'(x) (curves 1 and 2) 
electric field (curves 3 and 4) 
specific power (curves 5 and 6) 

In nitrogen, at atmospheric pres- 
sure, behind a foil of 50 y thick- 
ness for -electric, field Intensi- 
ties of 2 and 4 kv/cm, respec- 
. tlvely; Initial beam energy 200 p 
keV, beam current 5 • 10^ amp/cm'^. 
Dashed lines (7, 8) denote the 
levels of specific power corres- 
, ponding to the case of uniform 
Ionization. 


5 


Specific power generated in the gas discharge upon neglecting shock 
ionization is calculated from: ■ 


Pu^iE 


e|iC/o* 


) 'vr ' 


Xm«d. 


(13) 


(TVp'P(x) 

Equation (13) shows that specific power generated in the discharge in- 
creases with decreasing i.e., in the direction of . the anode. In the 
case of uniform ionization VCx) = specific power generated in the' 
discharge is equal to 


Po 

From (13) and (14) we obtain, . 

Pn 


ent/o*yy. 

(i*vr 


(14). 




If one assumes that H', 


P. Vy .'F(x) 

'V- 'i' (x = d), then ; 


(15) 


hence, we can 


In Figure 3 at cur- 


expect in the vicinity of the anode that 

4 2 “ . ® 

rent density 5 • 10 a/cm initial electron energy 200 ke'V, aluminum 

foil thickness 50 ym, U = 80 and l60 kv , respectively, we show the dis- 
tribution 'l'(x) curves 1, 2, E(x) - curves 3, Pjl - curves 5, 6. 
Dashed lines 7 and 8 show the level of specific power P^. For given 
Uq and d, shock ionization is negligibly small. From Figure .3, it is 
apparent that in the vicinity of the anode we not only get, an Increase 
in the field • - ^ but also, in the • generated specific 


power. Moreover, > P^j in spite of the fact that 
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. Tr» 

I Pn{x)dx , 

- — -4= < 1- 


(16) 


SMALL TIME NON-STATlONARY SOLUTIONS 


For short times 


1 




, and large Y we can neglect the drift 


and the recombination terms in (1) and in the discharge column we have 


The drift term is negligible if ' < '^'(x) for t = t , which 



Is equivalent to the Inequality 





dx 


( 17 ) 


where we set t = x . The system of equations (1) with the condition 
(17) gives an equation for the charge p: 


— p=eg-\-e\iEt—— 

dt t dx 


( 18 ) 


From the right-hand side of (l8), we see that the volume charge is 
created at the expense of the thermalized electrons :g and the gradient 
of the generation rate of the current carriers. In the regions where 
we can neglect the creation of volume charge at the expense/nf d^/dx, 
we obtain from ( 18 ):, 


Jexp[— af*(l— u)]du, O' 




(19.) 


2e 


The time dependence of the volume charge density p. shows an extre 
2 • 

mum... For at <<.1, the charge density grows in time as , 

2 




and for large t (but smaller x ), we obtain 

■ r 

eg 


2at 

where p is a decreasing function of time,. 


Let- us estimate themaximum p . The function ■ 

. Jexp[— cf*(l— u*)]du j 

0 ' ’ ■ ' 

is limited by the Inequality. , - - • - 

^exp (— af’)</<(l— exp(— />0. ; 

Hence, the maximum value 

p 

The left-hand side of the inequality has the extremum at at =0.5 
and the right-hand side at at^ 'v 1.25. Hence, :-.’0. 43 < max f < 
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— 1/2 

0.64. a is the extremum value of the density 


eg 

2Va. 


and the 


evaluation of the maximum field at the expense of thermalized elec- 
^trohs is 


9 J 


Pm(x) 


dx 


e . 


4 J 


Ud 


g{x)dx 

U(x) AR<^>e 


R>d.. 


( 20 ) 


= 1 a/cm^, e, = 35 ev, 3=2 
^ >8 
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10 cm'^/sec, we 


sec, Va = 5 • 10 sec 


experimental time a 


- 1/2 


For beam currents 
have: t = 5 • 10“' 

Sq 

= 2 • 10 ^ sec, and for distances of the order of R d, we obtain 
E,p 5 kv/cm. Because a it follows from (20) that Eip'v/J^. The 

above estimate shows that, under the conditions considered, the con- 
tribution of thermalized electrons to the field is appreciable for cur 

rents of the order of 1 amp and higher. Moreover, in a time propor- 
- 1/2 

tlonal to a the volume Charge density undergoes relaxation asymp- 
totically as t The ratio of the maximum value of charge density of 

thermalized electrons p to their saturation density is Indepen- : 

m n ■ 

dent of g and is determined by the species of gas under consideration ' 


"Pm' Tj/ep ' 
,P- ep* 


The relaxation is connected with the time growth of conductivity (if 

Let us consider next the building up of volume charge neglecting . 
thermalized electrons. Let us use the equation for. the total current 
which follows from (1) 


at z 6 


( 21 ) 


Assuming that the voltage along the discharge column remains prac- 
tically constant, after Integration with respect to, x, equation (21) 
gives ' . . ■ 


J{t)=-—\n{x,t)E{x,t)dx. 




_.:-:..:i22.) 


Given that at high electron pair generation rates 










( 23 ) 
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and assuming that y(x) is determined only by the average field—U^/d, 
and therefore is time-independent, from (.21), (22), and (23), we obtain 
an integral equation for the field: 


d ef^x„ 


.e-tJvvF'dx'je^'£(x'.V)th«y'F'M«'. (24) 


where i')=vlncht'V^f'p, v^en/e?, 


Equation (24) can be solved by- Iteration. As the zero-order ap- •; 

• i 

proximation, let us take E = U^/d. Then for t it follows from . ‘ 


(24) that 


I ywt}lty^f^dx 


Edx.t) ^ r I <yy> 

E - l-^m "v'F(z)’ 


(25) * 


As in equation (8), here E(x, «>) 't l/YXxJ. The results of (25) | 

agree with (8) bo within the constant coefficient, and are val.ld_.if i • 

-- r— / ^ 

* For* small times when n(x, t) = 'i'(x)t from (24), we | 
obtain --^ ■■ — i..,.- -4 , 


; E,{x,0 ; <'F> . / " 


'F(x) 




C26) !. 


l.e., for at >> 1, but t < Xpt 


Edx,t)._ <'F> . 

E ~ 'FUr 


(27) : 


.Comparison of (27) and (25) shows that the electric field intensity in 
the discharge column goes through an extremum. , , i 

As . is apparent from (25), the relaxation of the field from its ex- 

’.tremum to the stationary value is determined by recombination and takes 

place in a time t 'v. (2 - 3) .t . The field distribution at the time of 
/ , ■ P' ■ 2 ' ' ' 

thd extremum can be calculated exactly. For at •>> 1, > t, n(x, t), , 

I " dE' F 

■■h 'i'(x) • t, we can take — «0 and then it follows from (21) that 
/ , ■ . dt , ■ 


where c is a constant determined by the condition f Edx~Ut. 


Therefore, 


E / 1 




i<4> 


- 1/2 

This result is more exact, the better inequality t ^ >> a is satis- 
fied . . 


Specific power at extremum field (27) is determined by; 


p.. ^ V2epeP / 1 


when at >> 1, since 




(28) 


Thus, short-time considerations show that for non-uniform ioniza- 

- 1/2 

tlon after a time of the order of a , the field and the specific 
power undergo a significant redistribution within the gas volume. Af- 
terwards, during the title (2 - 3) t , one sees a decrease in inhomo- 
geneity. Since the first stage takes place in a relatively short’, time 
afterwards, conditions are especially favorable for those mechanisms 
of instability of volume gas discharge which are determined by the 
field intensity rather than specific power. ' , 

We must remark that (28.) does not exclude the possibility that 


'5. EXPERIMENTAL RESULTS - 


The fact that during electron injection the electric field in the'. 

gas gap can be strongly distorted is experimentally verified by the 

observation of instability of the gas discharge when the energy spec-. 

trum of the injected electrons is broadened [9, 10]. For injection 

-4 ■ ' 

pulse times of 10 sec., the. increase in the fall-off of. the voltage 
at the top of the accelerating voltage pulse always led with other 


10 - 


variables held constant to the Instability of the discharge [9]. The 
instability first appears in the/ anode region, i.e., specifically at 
the location where field. Increases must occur as a result of redis- 
tribution of the potential caused. by the inhomogeneous distribution of . 
conductivity in the gap.- We must remark, that the large duration of the 
pulse front of the accelerating voltage can likewise be a cause of in- 
stability. . 

The effect of the spectral characteristic of the electron beam of; 
duration 10“ sec. on the instability growth time was investigated in 
[10]. It was shown that lowering the maximum, energy of the injected 
electrons from 180 keV (current, 50 amps) to a value for. which the cur- , 
rent was. 1 amp, sharply reduced the instability growth time for dis- 
charge in nitrogen at 7 atmospheres. In a number of cases this reduc- 

_o 

tion was by a factor of IQO and meant a t_ 10 ^ sec.. .Such short in- ; 

n 

stability growth times are caused by streamer.-'discharge initiated in 
the region of field Increase, for which [11] ' ' ■ 

(29) 

av- ’ 

8 

where p is the gas pressure and 10 is the critical electron number 

-9 

in the avalanche. At nitrogen pressures of p = 7 atm and t 'v 10 , 

sec., the electrical field in the region of inhomogeneity growth must 
substantially exceed 10'^ volts/cmi Such a field can be easily due to , 

. thermalized, low-energy electrons [5]. We should note that formula 
(29) may be useful in the Investigation of fields due to electrons, 
thermalized in gases. 

In conclusion, the authors express their gratitude to Yu. D. Koro- 
lev .for helpful discussions. 
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